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We examine the role of space-time geometry in the non-adiabatic collapse of a star dissipating 
energy in the form of radial heat flow, studying its evolution under different initial conditions. The 
collapse of a star with interior comprising of a homogeneous perfect fluid is compared with that of 
a star filled with inhomogeneous imperfect fluid with anisotropic pressure. Both the configurations 
are spherically symmetric, however, in the latter case, the physical space t — constant of the config- 
urations is assumed to be inhomogeneous endowed with spheroidal or pseudo-spheroidal geometry. 
It is observed that as long as the collapse is shear-free, its evolution depends only on the mass and 
' size of the star at the onset of collapse. 
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I. INTRODUCTION 

Study of the evolution of a non-adiabatically collapsing spherical compact star under the influence of its self-gravity 



in the presence of dissipative forces and determining the final state it attains, is a field of active interest in theoretical 
^ ■ astrophysics. The Vaidya (ij metric appropriately describes the exterior space-time background of such collapsing 
^j^. stars. It has been generally expected that various factors contributing to the space-time geometry and matter-content 
I— I' within it, such as shear, inhomogeneity, local anisotropy, presence of electromagnetic field etc., may directly affect 
the evolution of the collapse. This is not surprising since the Einstein's field equations Gap = SttTq^ clearly imply a 
" ■ mutual correspondence between the composition of a gravitating body and geometry associated with its space-time. 

Banerjee et al [2] had presented a metric describing space-time of a spherically symmetric, shear-free, collapsing 
■ object in the presence of heat flux. Schafer and Goenner[3|] used this space-time metric and examined various features 
I of the evolution of the collapse of a shear-free homogeneous and isotropic configuration of matter. Pinheiro and 
Chan 0, @ analyzed the collapse of a matter distribution of uniform matter density with shear. We have examined 
the collapse of a relativistic star with matter density and fluid pressure decreasing radially outward from the centre, 
, considering the background space-time as having geometry different from that of a uniform density star. The space- 
rs) ' time is spherically symmetric and it's associated t = constant hyper surfaces have the geometry of either a 3-spheroid 
, or pseudo-spheroid. The use of space-time geometry of this type to obtain relativistic models of compact stars with 
^ ' no reliable information about their matter content was first proposed by Vaidya and Tikekar Q . Subsequently several 
• ^ ' physically viable models of realistic stars obtained adopting this approach have been discussed by many authors (e.g., 
^ " see Ref. 

In this paper, we have examined the impact of inhomogeneity of background space-time and anisotropic stresses on 
the collapse rate of a relativistic star which has lost its equilibrium. In Section |TT1 the field equations for a spherically 
symmetric anisotropic fiuid undergoing non-adiabatic radiative collapse on the background of a spherically symmetric 
inhomogeneous space-time have been formulated and the implications of prescribing the 3-spheroidal or 3-pseudo- 
spheroidal geometry for the interior physical 3-space of the collapsing configuration have been examined. In Section 
mil the nature of the space-time outside the collapsing star is discussed, the conditions to be fulfilled across the 
boundary surface separating it from the interior space-time have been formulated and an outline of the procedure for 
studying the evolution of the collapse is given. The concluding Section [TVl is devoted to discussion of our observations 
about the impact of the inhomogeneity on the evolution of the collapse of compact stars. 
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II. INTERIOR SPACE-TIME 



We consider a compact spherical star in equilibrium with interior space-time described by metric 

dsl = -A'oir)dt^ + Biir)dr^ + r^{d9'' + sin^ 0d<l)^), (1) 

wherein the coordinates have been labeled as x'^ — t, = r, x"^ — 6, x^ — (j). We assume that the interior content of 
the star is an imperfect fluid with energy-momentum tensor 

Tafi = [{p)s + {pt)s\UaUl3 + {ptjsOaP + [{Pr)s - {Pt)s]XaXP- (2) 

In Eq. ([2]), {p)s denotes the energy density, (pr)s and {pt)s respectively denote fluid pressures along the radial and 
transverse directions, u" is the four velocity of the fluid, while is a unit space-like four vector along the radial 
direction. The Einstein's field equations 

Ra/S - -^gapR = SnTaiS, (3) 

in system of units rendering G = c = 1, determine these dynamical parameters in terms of the metric coefficient 
functions as 

1 1 2A' 



r 



2 B^r^ rAoBl 



8TT{pt)s = 8TT[{pr)s+As{r)], (6) 

4" A' R' 4' R' 1 1 

^^"^^ A^Bl rAoBi vBl A^Bl Blr^^r^' 

In above As(r) represents the parameter which is a measure of anisotropy of the stresses. The subscript 's' is used to 
denote the dynamical variables associated with the compact star in equilibrium state. The interior space-time of the 
star after it begins to collapse on loosing its equilibrium will be described in general by the metric 

dsl = -A^{r,t)df + B^ir,t)dr^ + C^{r,t){d0^ + sin^ 0d(l)^). (8) 

The interior matter of the star will be accompanied with radial flow of heat as it shrinks in size and the energy- 
momentum tensor associated with it will be 

Tap = {P+ Pt)UaUi3 + Pt9al3 + {Pr " Pt)XaXlS + laUp + QpUa, (9) 

where is the heat flux vector. Einstein field equations ^ governing the evolution of the collapse is a highly non- 
linear system of second order partial differential equations. We adopt separability of metric variables which facilitates 
examining their implications. Without any loss of generality we couch the space-time metric of Eq. in the form 

dsl = -Al{r)dt^ + R^{t)Bl{r)dr'^ + r'^Y^ {t)(de'^ + sin^ edcj)"^), (10) 

and express the system of five field equations in the following form: 

1 , , 1 fy^ ry\ 1 / 1 1 \ 

Snp = ]^Mp). + ^(^ + 2— (11) 

8m = ir.MPt)s - T2 I §7 + ^ + - I ' (13) 

87rA(r,i) = -^B^A, + + (14) 

2 / A'Y R Y \ 
= + Z7^- 775^ ■ (15) 





AnR^Y rR^ tK^Y 
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In Eq. p4)) . the function A(r, t) represents evolution of the anisotropy in the collapsing configuration. The form of this 
function indicates that the evolution of anisotropy is linked with presence of shear (V ^ R) in the background. Even 
a compact star with a perfect fluid interior with vanishing (r) which starts shrinking in size when its equilibrium is 
lost will subsequently witness evolution of anisotropy in stresses as collapse proceeds. Apparently, shear is responsible 
for the development of anisotropic stresses in collapsing matter. 

The metric functions in Eq. ([T|) describing the interior space-time of a compact star in equilibrium are determined 
by solving the appropriate set of Einstein's field equations which needs reliable information about the equation of 
state of the matter content to be furnished a priori. For superdense stars with matter density exceeding nuclear 
density, this information is not available with certainty. In the absence of such information we adopt the alternative 
approach suggested by Vaidya and Tikekar (gj of prescribing a suitable geometry for the physical 3-space of the star 
and write the metric describing the interior space-time of the collapsing star as 



dst 



-Aoirydt^ + Rity 



1 



1 — kr^ 



dr^ +r^R^{t){de^ +1 



(16) 



where, we have set R{t) = Y{t). The geometry of the physical space t = constant hyper surface of the space-time 
described by (fT6)) is (i) spheroidal for k = -|-ve, and (ii) pseudo-spheroidal for k = — ve. Mukherjee et al [l^l and 
Tikekar and Jotania have shown that the respective matter distributions are physically plausible if A > 3/17 and 
A < — 1. The constant A is the measure of departure from homogeneous 3-spherical geometry of a uniform density 
star. When fc = 0, the physical 3-space is a flat Euclidean space. Substitution 



1 + Afcr2 

1 — kr'^ ' 

in Eq. ([7]) puts it in the form of a second order differential equation 



(17) 



(l + Afc-A/fcx^)^ 
dx'^ 



Xkx 



dAo 
dx 



Afc(Afc + l)- 



{l + \k- Xkx'^ySirAs 

fc(l-x2) 



(18) 



determining the metric parameter Ao{r) where, x denotes new independent variable defined as a;^ = 1 — kr^ . In terms 
of a new dependent variable 

^{x)^Ao{x){l + Xk-Xkx^)-^/'^, (19) 

Eq. IH]) assumes the form 

rf^^- r2Afc(Afc + l)(2AA: + l)-(4A/s + 7)A2fc2a;2 (1 + A/s - Afca;2)87rA, 



dx"^ 



* = 0. 



(20) 



4(1 + A/c - Afca;2)2 k{l - x^) 

^(a;) can be determined from Eq. (|20|) only when the nature of anisotropic parameter As is known in the interior 
of the star. It is expected that A^ should be regular everywhere in the interior. The equation (|20p admits a simple 
tractable regular solution 



= C + Dx, 

with C and D as arbitrary constants for 

fc(l - x^)[2Xk{Xk + l)(2Afc -I- 1) 



A., = 



i'iXk + 7)X'^Px'^ 



(21) 



(22) 



327r(l + Afc - Afca;2)3 

The above expression for As is regular at all interior points. The metric function Aoir) accordingly is explicitly found 
to be 



Ao{r) = (1 + Xk\^y^\C + DVl-kr^). 
Subsequently the interior metric for space-time of the shear-free collapsing star with matter will be 

'1 + Xkr^' 



ds^_ 



-(1 + A/cV2)i/2(c + Dy/l ~ kr^ydt^ + R^{t) 



1 — fcr2 



dr^ + R\t){de^ + i 



Space-time of Eq. ((24)) is an inhomogeneous generalization of the space-time of the metric of Banerjee et a/Q 

dr2 



dsl^ = -(C - Vl - kr'^fdt^ + i?2(t) 



1 — fcr2 



r^{de^ 



(23) 



(24) 



(25) 



as it can be recovered on setting A = and D = —\ from it. 

The progress of the collapse is directly linked with the metric function R{t) which plays the role of scale factor as 
the configuration shrinks in size. The determination of R{t) is discussed in the following Section. 
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III. EXTERIOR SPACE-TIME AND JUNCTION CONDITIONS 



Since the non-adiabatic collapse of the stellar configuration is accompanied with heat flux the space-time in its 
exterior is expected to witness presence of radiation and will be appropriately described by the Vaidya metric 



\d0' 



sm 



(26) 



with m{v) and E denoting respectively the total mass and the shrinking boundary of the star, both functions of the 
retarded time v. 

The metrics of the interior and the exterior space-time of the star should join across its boundary surface - a 
2-sphere of the evolving system. This is done by using the set of boundary conditions for a spherical collapsing matter 
distribution in the presence of heat flux completely formulated by Santos [l6j which ensure continuity of the first and 
second fundamental forms of the space-times. These conditions imply the following relations: 



= l- 



2m 



2^1 dv. 

dv , 



m(u) = 



[qR{t)Bo]^, 




(27) 
(28) 

(29) 

(30) 
(31) 



where, the functions Bo{r) and Ao{r) are as given in Eqs. ([TT)) and ([23)). respectively. The mass function ■m{r,t) 
defined as 



m(r, t) 



rR{t) 




(32) 



represents the total mass enclosed within the boundary surface r < at time t. It is expected that before the collapse 
sets in there will be no heat fiow across the boundary surface of the star. The corresponding boundary condition 
implies that {pr)s{'r = ''e) = 0. The condition pO|) leads to a differential equation 



2RR 



2nR = 0, 



(33) 



where, 



A' 



Bn 



(34) 



is a constant. Following Bonnor et al\Vi 



we write Eq. p3p as a first order differential equation 
R = — ^(1- V^), 



(35) 



which admits 



1 

t = - 

n 



R 

- + Vr + ln{l ~ y/R) 



(36) 



as its solution. It follows that at i — > — oo, the time when collapse is assumed to begin, R — 1 and i? = at < = 0. 
Hence R{t) will be a decreasing function of time. 

Let be the radius of the star when it was in equilibrium and be the corresponding mass of the star so that 
2ms/rs < 1. When instability sets in, the super-massive star loses its equilibrium and begins to collapse. We denote 
this initial epoch when collapse sets in as t = — oo so that R{t) = 1. The comoving boundary surface (ri?)x; = rsR(t) 
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TABLE I: Estimates of Rbh for initial stellar configurations of mass ms — 3.24 Mq and radius — 18 km with different 
compositions of interior stellar matter. 

^ C D k Rbh 

5 2.73837 -2.23899 0.00049 0.28196 

1.02725 -0.49999 0.00164 0.28196 

-5 -0.68817 -1.16291 0.00122 0.28196 



then starts shrinking until it reaches its Schwarzschild horizon value [rR{tBH)]Y: — fsRitsH) — 2m(w), with t — Ibh 
denoting the time of formation of the black hole. From Eq. ([29| . the mass of the evolving star at any instant will be 



m{v) = 



rrisR -i T2~(l 

^0 



RY 



where, 



(A + l)fcrg 
2(1 + Afcr2)- 



Consequently, the condition (tRbh)!! = 2m(w) yields 

Rbh = 



AnBo 



AnBo 



2ms 



and the time of black hole formation is obtained as 



1 

tBH = — 

n 



Rbh 



^R^+ln{l - ^JR BH) 



(37) 



(38) 



(39) 



(40) 



We note that Rbh depends on the initial values of the functions Aq and Bq on the boundary of the star when it starts 
moving inward as the collapse begins. The matter content of the star when the collapse begins is considered to be an 
inhomogeneous distribution of anisotropic matter which implies that A 7^ 0. If A = 0, < fc < 1, the star comprises of 
an incompressible distribution of perfect fluid with uniform matter density. The constants C, D, k and A appearing 
in the functions AqIt) and BqIt) are determined by the boundary conditions which join the interior space-time metric 
of the star with the Schwarzschild exterior metric describing its exterior space-time, at the boundary of separation 
. We stipulate the relevant boundary conditions as requiring continuity of metric tensor on Vs as 



1 



2ms 



1/2 



1 



2m, 



-1/2 



and vanishing of fluid pressure along radial direction (p,.)s(rs) — 0. These conditions determine 



C 



D^l-krli?, + A - Afc + 4Afc2r2 + X^k'^rl) 



(1 + A- A/c 



(41) 



(42) 



Note that Rbh depends on the values of Ao(r) and Bo{r) on the star boundary when collapse begins. The values 
of ^o(''^) and Bq^v) on the star boundary are completely determined by the mass and size of the star at the onset 
of collapse. Hence Rbh depends on the initial mass and size of the star. Therefore, in this framework, the time of 
formation of black holes Ibh and the radius vbh — TsRbh of the black hole remain the same for all stars of identical 
initial mass and size irrespective of the choice of A. We considered several stellar configurations in equilibrium of 
identical initial masses and radii (we assumed ms — 3.24 and radius Ts — 18 km) with different values for the 
parameter A implying different compositions of interior stellar matter and obtained estimates of Rbh using numerical 
procedures. We observed, as indicated in the Table HI that all configurations yielded almost the same Rbh- 

Accordingly, the parameter A seems to have no effect on the collapse rate. The inhomogeneous character of the 
interior matter content as implied by the geometry (spheroidal, pseudo-spheroidal or spherical) of interior physical 
space and nature of matter composition (fluid with isotropic or anisotropic stress structures) have no impact on the 
collapse rate of a stellar configuration. The function R{t) is only a scaling parameter and carries no information about 
the geometry or material composition of the star in this set up. This is consistent with implications of Birkhoff 's 
theorem according to which no observer located in the exterior region of a spherical star can get any information 
about the interior state of the star except its mass and size. 
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IV. DISCUSSION 

Considering the non-adiabatic collapse of a spherical star of uniform density Oliveira et al^^ have shown that 
some of the features like time required to form the black hole depend on the initial mass and size of the star. In our 
work, we have generalized the space-time metric of Banerjee et al to investigate the effects of inhomogeneity and 
anisotropy on the collapse. Our study indicates that collapse rate for stars having identical initial mass and size does 
not depend on the geometry of their interior physical space. They take same time to collapse into their horizon and 
become black holes. Since the geometry of the physical space of a star in equilibrium is related with the composition 
of the interior matter of the star, this implies that stellar configurations having identical initial masses and radii 
irrespective of their composition, after loss of equilibrium, will take same time to become black hole. 

The analysis above is based on the following two assumptions. The space-time in the interior of the star (1) is 
spherically symmetric and shear-free, and (2) described by metric with metric functions as separable functions of space 
and time variables. Analysis based on such simplifying assumptions may give the same result unless new techniques 
are invoked to extract information about the change in the equation of state of the interior matter of a collapsing 
relativistic star and the mechanism for energy transport. The mass and radius of a spherical star in equilibrium will 
be the only governing parameters which determine its collapse rate and time for its becoming a black hole after it 
begins to collapse unless shear develops in its space-time structure. Therefore, it is essential to consider non-adiabatic 
collapse of stellar configurations on the background of spherically symmetric space-time admitting shear R{t) ^ Y{t) 
and study the impact of initial anisotropy and inhomogeneity on the evolution of the collapse. 
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